Computing absolute free energies of disordered structures by molecular simulation

We present a Monte Carlo simulation technique by which the free energy of disordered systems can be computed directly. It is based on thermodynamic integration. The central idea is to construct an analytically solvable reference system from a configuration which is representative for the state of interest. The method can be applied to lattice models (e.g., the Ising model) as well as off-lattice molecular models. We focus mainly on the more challenging off-lattice case. We propose a Monte Carlo algorithm, by which the thermodynamic integration path can be sampled efficiently. At the examples of the hard sphere liquid and a hard disk solid with a defect we discuss several properties of the approach

A density-functional theory study of the confined soft ellipsoid fluid

A system of soft ellipsoid molecules confined between two planar walls is studied using classical density-functional theory. Both the isotropic and nematic phases are considered. The excess free energy is evaluated using two different Ansätze and the intermolecular interaction is incorporated using two different direct correlation functions (DCF’s). The first is a numerical DCF obtained from simulations of bulk soft ellipsoid fluids and the second is taken from the Parsons–Lee theory. In both the isotropic and nematic phases the numerical DCF gives density and order parameter profiles in reasonable agreement with simulation. The Parsons–Lee DCF also gives reasonable agreement in the isotropic phase but poor agreement in the nematic phase